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ABSTRACT 

Black hole mass (Mbh) is a fundamental property of active galactic nuclei (AGNs). In the distant 
universe, Mbh is commonly estimated using the Mgll, H/3, or Ha emission line widths and the 
optical/UV continuum or line luminosities, as proxies for the characteristic velocity and size of the 
broad-line region. Although they all have a common calibration in the local universe, a number of 
different recipes are currently used in the literature. It is important to verify the relative accuracy and 
consistency of the recipes, as systematic changes could mimic evolutionary trends when comparing 
various samples. At z = 0.36, all three lines can be observed at optical wavelengths, providing a 
unique opportunity to compare different empirical recipes. We use spectra from the Keck Telescope 
and the Sloan Digital Sky Survey to compare Mbh estimators for a sample of nineteen AGNs at this 
redshift. We compare popular recipes available from the literature, finding that Mbh estimates can 
differ up to 0.38 ± 0.05 dex in the mean (or 0.13 ± 0.05 dex, if the same virial coefficient is adopted). 
Finally, we provide a set of 30 internally self consistent recipes for determining Mbh from a variety of 
observables. The intrinsic scatter between cross-calibrated recipes is in the range 0.1 — 0.3 dex. This 
should be considered as a lower limit to the uncertainty of the Mbh estimators. 

Subject headings: black hole physics: accretion — galaxies: active — galaxies: evolution — quasars: 
general 



1. INTRODUCTION 

Understanding the growth of supermassive black holes 
along with their host galaxies is one of the ftindamental 
quest i ons in current astr ophysics (e.g. iDi Matteo et al.l 
[2005t iCroton et all 120061) . Black hole mass (Mbh) is 
a key parameter in revealing the nature of black hole- 
galaxy coevolution as well as the physics of active galactic 
nuclei (AGNs). However, direct mass measurements us- 
ing the motions of gas and stars in the sphere of influence 
of a cen tral black hole is limited to very nearby galax- 
ies (e .g. iKormendv fc Gebhardtl [2OOII : iFerrarese &: FordI 
[2OO5I) . 

Beyond the very local universe, the so-called "virial" or 
"empirically calibrated photo-ionization" method based 
on the reverbera tion sample is popularly used for ac- 
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Wandel et ai]|1999t iKaspi et al.ll2000l 
|2006( ). This method utilizes broad 



line widths as velocity indicators and monochromatic 
continuum or line luminosities as indicators of broad- 
line region size, hence estimating virial Mbh- A com- 
bination of the Mgll, H/3, or Ha broad emission line 
widths and the 3000A, 5100A, H/3, or Ha luminosi- 
ties is typically used, depending on the redshift of 
the source and the observational setup. Several equa- 
tions have been presented in the literature to esti- 
mate Mb h using va rious combinations of these indica- 
tors (e.g., pVoo fc U rrv 2002a, b; McLure lTjarvigl 120021: 
iTreu et all 120041 : [KoUmeier et al.l 120061 : iGreene fc Hoi 
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Although all three emission lines have a common cal- 
ibration based on the reverberation sample in the local 
universe, it is important to verify that different recipes 
give consistent results; any systematic changes could 
mimic evolutionary trends given that different recipes are 
often used in various studies. 

At z = 0.36, all three lines can be observed at optical 
wavelengths, providing a unique opportunity to cross- 
calibrate the different methods of Mbh estimation. Using 
data from the Keck Telescope and the Sloan Digital Sky 
Survey for a sample of nineteen AGNs at z = 0.36, we 
compare the different methods of estimating Mbh, and 
derive a set of self-consistent equations for Mbh estimates 
using every combination of velocity scale (FWHM and 
line dispersion (June of Mgll, H/3, or Ha) and luminosity 
(3OOOA, 5IOOA- nuclear and total - H^, or Ha). 

The paper is organized as follows. In § [2] we describe 
the sample selection, observations, and data reduction. 
In § [3] we describe our line fitting process, based on ex- 
pansion in Gauss-Hermite series, and the resulting lumi- 
nosity and width measurements. In § [4| we review the 
various formulae adopted in the literature, and compare 
the various Mbh estimators. In §[5] we present our self- 
consistent recipes. Section [5| summarizes our results. 

Throughout this paper magnitudes are given in the AB 
scale. We assume a concordance cosmology with matter 
and dark energy density f2„i — 0.3, f2A — 0.7, and Hubble 
constant Ho=70 kms~^Mpc^^. 

2. SAMPLE SELECTION, OBSERVATIONS AND DATA 
REDUCTION 

The AGN sample was initially selected for stellar veloc- 
ity dispersion (cr* ) measurements to measure the Mbh-o'* 
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Fig. 1. — Flux-calibrated spectra. The SDSS spectra are shown in black, and the Keck spectra are shown in blue and red. The Mgll 
line can be seen on the far left of each wavelength range, while H/3 is located in the center and Ha to the far right. 



relation at z = 0.36 ()Treu et al.ll200l I Woo et al.ll2006f l. 
Readers are referred to the papers by Woo et al. (2006) 
and Treu et al. (2007) - where Keck red spectra (5900A- 
7500A) and Hubble images of the sample were presented 
- for more details. The relevant properties of the ob- 
served objects are listed in Table [TJ 

High signal-to-noise ratio spectra of nineteen targets 
were obta ined with the L ow Resolution Imaging Spec- 
trometer ()Oke et al.lll995l hereafter LRIS) at the Keck- 
I telescope in five r uns between Marc h 2003 and July 
2005, as detail ed bvlWoo et all (I2006D . The red setup 



is described bv lWoo et all (|2006f ). In the blue, the 600 



culating synthetic g' and r' magnitudes from the spectra 
taking into account the SDSS bandpasses and finding the 
constant multiplicative factor appropriate to match the 
SDSS photometry. The final flux-calibrated spectra of all 
nineteen observed AGNs are shown in Figures [T] and [2j 

3. MEASUREMENTS 

The relevant quantities for our purpose are line widths 
and line and continuum fluxes. In this paper, line widths 
are measured as FWHM or line dispersion, i.e. the square 
root of the second central moment: 



lines mm grism was used, yielding a pixel scale of 
0.63AxO'.'135 and a resolution of ~145 km s'^ Note 
that objects S16, S3 1, and S 99 included in the papers by 
IWoo et al.l (|2006l l or lTreu et a l. (2007) lack Keck and/or 
SDSS spectra and are therefore not considered here. 

The reduction of the blue spectra was very similar t o 
that of the red spectra described by IWoo et aD (|2006f l , 
except that arc lamp emission lines for Hg and Cd were 
used for wavelength calibration due to the paucity of sky 
lines. The flux was calibrated using spectrophotometric 
stars or AOV type Hipparcos stars. Galactic extinction 
correction was applied to all data based on th e average 
extinction law derived in ICardelli et alj (|1989f) . The fi- 
nal step in the reduction process involved normalizing 
all spectra to the proper AB magnitude values from the 
Sloan photometric database. This was achieved by cal- 
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/(FA-C)(A-Ao)^dA 
/(Fa - C)dA 



(1) 



where F;^ is the flux density, C is the continuum, A 
is the wavelength and Aq is the central wavelength of 
the line. As described in this section, we derive line 
widths and fluxes by fitting a Gauss-Hermite series to the 
data after subtraction of the continuum, Fe emission, and 
when necessary a narrow component of the line. While 
the Gauss-Hermite fit is not necessary for the high S/N 
ratio Keck data - whe re the quantities can be measured 
directly from the data (|Treu et al.ll2004l : IWoo erani2006l ) 
- the fitting procedure is needed to obtain robust mea- 
surements for the noisier SDSS data. As we demonstrate 
below by comparing the H/3 measurements from Keck 
and SDSS data, this fitting procedure gives consistent 
results between the two datasets, and therefore indicates 
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Fig. 2.— As in Figure[T]for objects Sll to S29. 
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Fig. 3. — Example of nuclear Fe subtraction. The top panel 
shows a typical observed spectrum (black histogram) , together with 
the nuclear iron emission template matched in intensity and res- 
olution (red line). The bottom line shows the residuals after Fe 
subtraction. Note that the bump redwards of Mgll emission has 
disappeared in the residuals. 



that the quality of the SDSS data is adequate to measure 
the width and flux of H/3 and, especially, Ha, since the 
latter is considerably stronger. 



3.1. Gauss-Hermite Fitting 

The emission lines, especially for H/3, are often asym- 
metrical, thus making a symmetrical Gaussian ap- 
proximation of the line profiles undesirable. To ac- 
count for the asymmetries in the emission lines, we 
fit a truncated Gauss-Herm ite series to the profiles 
(Ivan der Marel fc Franxl [l993l ) . The main advantage of 
the Gauss-Hermite expansion is that it provides an or- 
thonormal basis set and that the coefficients of the Her- 
mite polynomials (commonly referred to as h^, /14, etc.) 
can be derived by straightforward linear minimization, 
leaving only two non-linear parameters (the center and 
the width of the Gaussian). Furthermore, the coefficients 
can be interpreted i n terms of the k inematics of the trac- 
ing population (e.g. lGerhardlll993[) . The best fit profiles 
arc then used to measure the luminosity, FWHM, and 
(Jiine of the emission lines as parameters for the Mbh 
formulae. 

We begin the fitting process with continuum subtrac- 
tion. We identify the continuum level on each side of 
the three emission lines, using a window of 60 A. In the 
case of MgH, a narrower 40 A window is used because 
the blue continuum of Mgll is close to the end of the 
observed spectral range. For MgH, the ranges typically 
used are 2660-2700 A for the blue continuum level and 
2930-2970 A for the red continuum level. Similarly for 
H^, the ranges are 4670-4730 A and 5080-5140 A for 
blue and red, respectively, and for Ha the ranges are 



4 



McGill et al. 




Fig. 4. — Gauss-Hermite broad line fits. From left to right, the 
boxes show the fits for the Keck Mgll line, Keck H/3 line, SDSS H/3 
line, and SDSS Ha line for each object. The continuum-subtracted 
line is shown in black, the broad component is shown in blue, and 
the broad-line fit is shown in red. 
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Fig. 5.— As in FigureHfor objects S07 to S12. 

6290-6350 A and 6700-6760 A. We then independently 
subtract the continuum by linear interpolation for each 
line. 

Together with the featureless continuum, we also re- 
move broad nuclear Fe emission underneath Mgll and 
up. For this task we use template spectra of I Zw 
1 kindly provided by Todd Boroson and Ro ss McLure. 
The p rocedure is similar to that described bv lWoo et al.l 
(|2006f ). A typical example of Fe subtraction in the wave- 
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Fig. 6.— As in Figure|4]for objects S21 to S29. 



length region around Mgll is shown in Figure [31 Re- 
moving nuclear Fe emission changes the measured width 
of H/3 by a negligible amount (the FWHM is unchanged 
and (Tiino is reduced by 0.013 dex), but it has a signifi- 
cant effect on Mgil (the FWHM is reduced by 0.027 dex, 
while iJiine is reduced by 0.106 dex). 

Before fitting the broad lines, we subtract out the 
narrow lines, ex tending the procedure described by 
IWoo et ahl (|2006f) . For H/3 this involves subtraction of 
the [OHI] A4959 and A5007 narrow lines, as well as the 
narrow H/3 line. We subtract [OHI] A5007 directly, and 
subtract [OHI] A4959 by dividing [OIII] A5007 by 3 and 
blueshifting. The narrow component of H/3 is subtracted 
by rescaling and blueshifting [OIII] A5007. The line ra- 
tio H/3narrow / [O III]A5007 was allowed to range between 
1/20 and the maximum value consistent with the absence 
of "dips" in the broad component (typically 1/10-1/7; 
e.g. Marziani et al. 2003). The adopted value of the scale 
factor is listed in Table [2] The narrow component of Ha 
is subtracted by multiplying the determined H/3 narrow 
component by 3.1 (Malkan. , 1983.; Oster brock .1989) and 
redshifting. 

No attempt is made to remove narrow Nil emission 
lines around Ha, since they are effectively rejected by 
the Gauss-Hermite fitting procedure as noise spikes, nor 
the narrow component (if present) of Mgll. 

The final step in approximating the line profiles in- 
volves fitting the Gauss-Hermite series to the broad lines. 
The fitting procedure finds the minimum by increas- 
ing the order of the Hermite polynomials as required by 
the data (i.e. only if the goodness of fit, measured by 
the reduced x^, improves). Most objects required up to 
order 6 polynomials (i.e. he) to plateau in reduced x^- 
Figures[4]to[6]show our best fits to the Mgll, H/3, and Ha 
lines for the nineteen objects. The derived measurements 
of (Tiino and FWHM, after removal of the instrumental 
resolution, are listed in Tabled 
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3.2. Luminosities 

The formulae for Mbh estimates require line luminosi- 
ties or monochromatic continuum l uminosities a t given 
wavelengths. For consistency with , Greene fc Hoi ()2005f ) 
the line luminosities are calculated from the total flux for 
the combined broad and narrow components of the H/3 
and Ha emission lines, where the broad component is 
taken as the best Gauss-Hermite fit, and the narrow com- 
ponent is taken as the appropriately scaled and shifted 
[OIII] A5007 narrow line (see § [SH]). We note that the 
narrow components contribute only a small fraction of 
the total flux of the Balmer lines (e.g., ~5% for H/3 - see 
Table 2, and similarly to Ha) and therefore they make 
a contribution of about 0.01 dex to the broad-line size 
estimates. The Ha to H/3 flux ratio ranges between 3 
and 7 as expected for Seyfert Is fe.g. lLacv et al.|[r982D . 

The total continuum luminosity at 3000 (5100) A is 
calculated from the average flux in the 2950-3050 (5050- 
5150) A rest frame. In this paper, we use the term total 
continuum luminosity to indicate the total luminosity as 
measured within the spectroscopic aperture, i.e. without 
removing the host galaxy contamination. The fraction of 
host galaxy contamination depends on the properties of 
the individual object as well as on the instrumental setup, 
and it is hence a source of scatter. Nevertheless, the total 
luminosity is often the only measurement available and 
therefore it is important to investigate estimators based 
on this quantity. 

The total continuum luminosities at 5100 A measured 
from the SDSS spectra agree to within a few per cent of 
those inferred fro m the Keck spectra and those listed in 
the paper by Woo et all ()200 6^. By comparing the val- 
ues lis ted here with respect to those given bv lWoo et al.l 
(|2006f ). we infer 0.014 dex as the error on the contin- 
uum (L3000 and L5100). For line luminosities, we com- 
pare measurements on the fit with measurements on the 
data, and we take the r.m.s. scatter as the average error. 
This error is 0.011 and 0.062 dex on Lh,3 and Lho, re- 
spectively. Nuclear luminos ities are taken fro m the HST 
measurements presented bv lTreu et al.l ()2007[ ). except for 
three objects (Sll, S28, and S29), where the nuclear lu- 
minosity has been estimated from scaling the total lu- 
minosity by the average nuclear fraction for the sample, 
0.31. In addition to measurement errors, AGN variabil- 
ity effectively limits the accuracy of the calibration of 
luminosity-based esti mators, to the typical le vel of vari- 
abilit y of 5-10% (e.g., iWebb fc Malkan 2000: .Woo et all 
l200l . 

3.3. Line Widths 

The Mbh estimators depend either on CTune or on the 
FWHM as a measurement of line width, and hence of 
the broad line kinematics. Both of these quantities were 
measured from the Gauss-Hermite fit to the broad lines. 
Errors on the Keck line width measure ments are obtaine d 
by comparison with those reported bv lWoo et al.l (|2006D . 
which were measured independently and did not rely on 
Gauss-Hermite expansion. The average error is 0.017 
dex for the Mgll and H/3 FWHM and anno- The errors 
on the SDSS Ha line width measurements are obtained 
by comparing the results for H/3 with the Keck fit and 
assuming that the error is the same for Ha and H/3. The 
average error is 0.051 dex for au^c and 0.040 dex for the 
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Fig. 7.— Distribution of FWHM to aunc ratios for Mgll, H/3, 
and Ha. For comparison, the expected value for a Gaussian is 2.35. 
The typical errors on the measurements are 0.03 dex for MgH and 
H/3 and 0.06 dex for Ha. 

FWHM. 

The distribution of the FWHM to cruno ratios is shown 
in Figure [T] The average ratio for H/3 and Ha is close 
to the Gaussian value (2.35), alt hough with larg e scat- 
ter, consistent with the sample of lPeterson et al.r (|2004). 
For Mgll the average is considerably smaller indicating 
significant departure from Gaussianity. A large range 
of FWHM/cTiine ratios indicates a large scatter between 
Mbh based on FWHM and Mbh based on aune for 
the sample since velocity is simply derived either from 
FWHM or tJiinc by multiplying by a constant as sh own 
in §|4](see detailed discussion bv I Collin et al.l (I2006D V 

Figure M compares the FWHM of Mgll with that 
of H/3 and Ha. The FWHM are correlated albeit 
with substantial scatter. Comparing all the veloc- 
ity scales, the average ratios and r.m.s. scatters 
(in parenthesis) are: (log(FWHMM.ii/FWHMH/3)) = 
0.02 ±0.03(0.13), (log(FWHMH/3/FWHMHc)) = 0.09 ± 
0.02(0.07), (log((7Mgii/(TH;3)) = 0.13 ± 0.02(0.10), 
(log((7H/5/crH«)) = 0.07 ±0.02(0.10). 

Summarizing these relations, the width of Ha is gen- 
erally narrower by ^2 0% than tha t of H/3 as expecte d 
from other studies fe.g. lShudeij[l984HGreene'fc Hal2005f) . 
while Mgll and H/3 are similar in FWHM but not in 
(Jiinc- These differences are expected given that differ- 
ent lines trace different species and that their shapes re- 
flect the ionization and excitation variations throughout 
the region. This finding implies that each line has to 
be calibrated independently as a velocity estimator and 
that one cannot go from FWHM to criinc using the sim- 
ple scaling for a Gaussian distribution. Therefore, we 
conclude that in general Mg II and Balmer line width 
cannot be used interchangeably, although in the case of 
the FWHM of Mg and H/3 the ratio is close to unity (see 
iMcLure fc Jarvij2002l : [Salviander et al.ll2007l for further 
discussion). 

As far as Mbh estimators are concerned, the scatter is 
of order 0.1 dex (i.e. significantly larger than the mea- 
surement errors), which sets a lower limit of ~ 0.1 — 0.2 
dex on the relative uncertainty of the cross calibration of 
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Fig. 8. — Comparison of the width of Mgll with that of H/3 
and that of Hq. The average errors are presented in the bottom 
right Conner. The two objects furthest away from the ratio=l hne 
(dashed line) are labeled for easy comparison with the line profiles 
shown in Figures [5] and |6] 

4. REVIEW OF OPTICAL-UV Mbh ESTIMATORS 



simple velocity estimators based on line widths. 

We begin this section with a list of the 12 formulae 
for Mbh estimation considered in this paper (§ 14. ip . In 
these formulae, we adopt a notation where Lsioo.t =total 
luminosity \L\ at A = 5100A, \jzwo,n ^nuclear luminos- 
ity \Lx at A = 5100A, and L3000 = AXa at A = 3000A. 
All formulae are given in the original notation, without 
applying any correction for different assumptions on the 
virial coefficient. 

In § l4.2l we will compare the various estimators to infer 
how much they differ when applied to the same sample 
of objects before presenting our cross-calibrated recipes 
in §13 

4.1. Summary of existing recipes 

iMcLure k Jarvid (l2002f ). iKollmeier e\~^ (|2006( l and 
ISalviander et alJ (|2007| ) give equations based on the 
width of the Mgll line and the optical/UV continuum 
luminosity: 



Ml 



M 



3.37 



L 



3000 



(i037 W 



0.47 



/ FWHMMgii 
V km s~^ 



M, 



0: 



(2) 



Mk = 2.04 ( .V" PWHM-^"Vm„ 



1044 ergs- V V kms^i 



(3) 



Msa = 10^ 



£5ioo,t / FWHMMgn y ^ 

1044 erg s-V V3000 kms-V ®' 



iGreene fc Hoi |2005) and lVestergaard fc PetersonI ()2006f ) 
present equations based on the width and luminosities of 

the H/3 and Ha broad lines: 



M, 



3.6x10'' 



1042 erg s-i 



0.56 



/ FWHMh;3 ' 

UOOO km s"\ 



Mo, 
(5) 



Mgo = 2.0x10^ 



1042 erg s-i 



0.55 



/ FWHMho 
V 1000 km s" 



M0, 



Mv/3 = 10' 



,6.67 



Shields et al.l 



1042 erg s-i 
(12001 . 



/ FWHMh/3 ^ ^ 
UOOO km s"^ 



(6) 
(7) 



iCreene fc Hoi ([2005) 
Vestergaard & Peterson' (2 0m . |Woo et al.' (2006) 



Netzer & Trakhtenbrot (200j7 and Trcu et al. (200?!) 
adopt the following formulae based on L5100 and the 
width of the H/3 broad line: 



M; 



Sh 



10 



7.69 



L 



5100, t 



1044 erg s- 



0.5 



/ FWHMh/3 
V3OOO km s"^ 



Mo, 
(8) 



Mg51 = 4.4x10^ 



Mv = 10'' 



^5100, n 



V1044 



erg s 



0.64 



FWHMh/3 
1000 km s"^ 



Mo, 
(9) 



L 



5100, t 



1044 erg s-i 



/ FWHMh/3 ^ ^ 
V 1000 km s" 



M 



Mw = 2.15x10** 



AfN = 1.05x10** 



^5100, t 



1044 erg s-i 



5100, t 



V 1046 erg s- 



0.69 



0.65 



Mt 



10' 



S.58 



L 



5100, n 



1044 erg s-1 



0.518 



3000 km s" 



FWHMh/3 
1000 km s"^ 



3000 km s" 



0) 
(10) 

Mo, 
(11) 

Mo, 
(12) 



Mo. 
(13) 



Note that the formulae listed above adopt different es- 
timators of broad-line region velocity (uiine or FWHM) 
and size (continuum luminosity or line luminosity). 

4.2. Comparison of Mbh estimators 

In this section we assess the relative consistency of the 
Mbh estimators taken from the literature by comparison 
with our fiducial black hole mass Mt (as listed in the 
paper by Treu et al. 2007). The choice of the fiducial 
estimator is based on reverberation mapping studies of 
local AGNs, which are mostly based on H/3 and £5100- 
These stu dies show that aiine is the most robust v elocity 
estimator (iPeterson et 311120041 : ICollin et al.ll2006D . They 
determine the slope of the size-luminosity rel ation, taking 
into account the host galaxy contamination (|Bentz et al.l 
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difference (logAfan — log Mt) and r.m.s. scatter in dex are listed in the bottom right corner. The formulae are taken directly from the 
literature with no adjustments for the difference in the choice of virial coefficients. 



l2006f ). and they set the virial coefficient by requiring 
that the Mbh-c* relation be the same fo r active and 
quiescent galaxies , modulo selection effects (jOnken et al.l 
[200l iLauer "eralll2007) . 

It is important to notice that the formulae Mqsi, Mq/j, 
Mca, Msh, and Msa have been calibrated on the isotropic 
spherical virial coefhcient (f=3/4 in the notation of Net- 
zer, 1990), Mm and Mk adopt f=l, and Mt, Mw, My, 
My^, and Mn are based on the recalibration of the virial 
coefficient given bv lOnken et al.l (|2004( )^. Therefore, we 
expect the first set of Mbh estimators to give lower val- 
ues than Mt by log 1.8 = 0.255 dex, and Mm and Mk 
to give lower values by log(1.8 x 3/4) = 0.130 dex. In the 
case of Msh and Msa, an isotropic spherical virial coeffi- 
cient was used. However, their equations were based on 
a different size-luminosity relation, hence making these 
formulae approximately equivalent to having the same 
virial coefRcient as Mt, Mwi and Mn (jSalviander et al.l 
I2007f) . As noted bv iTreu et all (|2007[ ). Mw, Mn, Msh, 
and Msa agree on average to within a few per cent of 
Mt. The more discrepant black hole masses are those 
with different virial coefficients, which can differ on av- 
erage by as much as 0.38 ± 0.05 dex (i.e. more than a 
factor of two) . Even renormalizing these formulae to the 

^ It is also generally assumed that FWHM=2(Tiine- This is con- 
sistent on average with our result for the Balmer lines, albeit with 
large scatter, but not for the Mg II line. 



same virial coefficient would still leave discrepancies of 
order 0.1 dex: after renormalization. My and Mq/s are 
approximately 0.1 dex larger than Mt while Mqsi, Mm, 
and Mk are still 0.1 dex smaller than Mt. This is ap- 
proximately twice the expected error on the mean given 
the size of the sample. 

These results show that systematic errors as large as 
0.38 ±0.05 dex can be introduced when comparing Mbh 
estimates based on different diagnostics or when com- 
paring AGN Mbh estimates to local samples with direct 
Mbh measurements from stellar or gaseous kinematics. 

We conclude by discussing the absolute c alibration of 
our fi ducial mass estimator. A study by ICollin et al.l 
()2006l ) finds that for the local sample with reverberation 
based Mbh , the widths of the Balmer lines measured on 
mean spectra are systematically larger than those mea- 
sured on the r.m.s. spectra for variable objects, suggest- 
in g that a smal l er vir ial coefficient than that advocated 
bv lOnken et"al] (|2004[ ) and used for our fiducial mass es- 
timator, should be adopted when measuring widths from 
mean spectra. This could possibly indicate that our fidu- 
cial Mbh are overestimated by approximately 0.15 dex. 
To investigate this effect, in Section [5] we apply our cal- 
ibrated recipes - based on our fiducial estimator - to 
the local sample with reverberation based Mbh, using 
published line widths and fluxes measured from single 
epoch and mean spectra. As discussed in the next sec- 
tion, we find an excellent agreement with the mass in- 
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ferred from reverberation mapping based on line widths 
from the r.m.s. spectra, indicating that no such cor- 
rection is necessary. Data for a larger number of objects 
with reverberation-based masses are needed to determine 
the zero point of the virial scalings more accurately, as 
discussed in the next Section. 

5. A SET OF SELF-CONSISTENT RECIPES 

In this Section we examine all possible combinations 
of velocity and flux estimators to produce a set of cross 
calibrated recipes. By comparing with our fiducial mass 
estimator Mt, we compute the r.m.s. scatter of the resid- 
uals to infer a lower limit on the intrinsic uncertainty of 
each recipe. 

In practice, we adopt the following relation: 

logMBH = a + 21ogviooo +/31ogL, (14) 

where viooo is a velocity estimator in units of 1000 km 
s""'^, L is a luminosity estimator in units of 10'*'' erg s^^ 
or 10^^ erg s~^, respectively for continuum or line lu- 
minosity, and we find the a. that best matches the Mt 
fiducial estimates. Our range in luminosities is too small 
to fit for /3 as well, and therefore, we assert the follow- 
ing fiducial values: /3 = 0.47 for L3000, 0.518 for Lsioo,™, 
0.67 for Lsioo,*, 0.56 for H/3, and 0.55 for Ha. These 
choices are based on the most current calibration of the 
size-luminosity relation for each wavelength / line. In par- 
ticula r, following the results of the study bv lBentz et al.l 
()2006l ) we adopt 0.518 as the slope of the broad-line re- 
gion size vs nuclear luminosity relation, while for the 
size vs total luminosity (i.e. including host galaxy con- 
tamination within the spectr oscopic aperture) r elation 
we adopt the slope given by iKasDi et all ^OM)- Al- 
though the former slope is to be preferred when the nu- 
clear luminosity is available, we also provide results for 
the second slope, which appears to be the best estimate 
whenever the light from the nucleus and from the host 
galaxy cannot be disentangled. In general, extrapola- 
tions well outside the range considered here are to be 
done with caution, since most of the local calibrators for 
the size-luminosity relation are Seyferts and PG quasars 
in 10^2 < L5ioo(ergs-i) < 10*6. 

We emphasize that this procedure produces self consis- 
tent mass estimates, but these all share a common uncer- 
tainty in the zero point. In practice, all a can be shifted 
by a constant if it turns out that a different value of 
the virial coefficient for the local sample of reverberation 
mapped AGNs is to be preferred. As a sanity check, we 
applied our calibrated recipes to the local sample with re- 
verberation Mbh (Peterson et al. (2004) and recent up- 
dates bv iBentz et al.l ([20061), iDennev et^ (|2006[ ). and 
iBentz et al.l (|2007D ). Based on the single epoch m ea- 
surements given by IVestergaard &: PetersonI ()2006D we 
compared our two estimators based on the FWHM of 
II/3, and on Lsioo,* and H/3 line flux, with the reverber- 
ation masses. The agreement is excellent, with average 
A log Mbh.icv - A log Mbh equal to 0.009 dex and 0.023, 
respectively, with r.m.s. scatter of 0.46 dex and 0.49 dex. 
Ba sed on the measur ements from mean spectra given 
bv ICollin et al.l ()2006l ) we compared our two estimators 
based on the FWHM and line dispersion of H/3, and on 
LsioOjtj with the reverberation masses. The agreement is 
again excellent, with average A log Mbh, rev — A log Mbh 



equal to 0.004 dex and 0.018, respectively, with r.m.s. 
scatter of 0.35 and 0.29 dex. The scatter of the difference 
is combination of effects from uncertainties in time-lag, 
line width, and luminosity measurements, and scatter in 
the size-luminosity relation. 

The best fit values of a together with the r.m.s. scat- 
ter are listed in Table [31 The best fit relations are shown 
in Figure [TOl We note that we included in the recipes 
the combination of H/3 criinf, and LRi nn.n. that was used 
to compute Mt bv lTreu et all ([20071). The goal of this 
exercise is to estimate the measurement errors associ- 
ated with the relation, given that the input parameters 
were measured independently for this paper, based on 
the Gauss-Hermite polynomial expansion fit. Since mea- 
surement errors such as narrow line and continuum sub- 
traction dominate over pixel noise, given the high S/N 
of the Keck data, the r.m.s. of these residuals divided by 
V2 is effectively the total measurement error, i.e. 0.03 
dex. Thus we can conclude that, for all practical pur- 
poses, the r.m.s. scatter that we observe for the other 
recipes is intrinsic scatter in the relation and not mea- 
surement error. This is also the rationale for not showing 
measurement error bars in the plots. 

Looking at Table [3l it appears that the smallest rela- 
tive scatter is obtained when comparing Mt to black hole 
masses based on the same velocity scale (Th/3 . This is not 
surprising, as the velocity scale enters with the square in 
the Mbh estimate and we have seen that velocity scales 
have typical relative scatter of 0.1 dex. However, the 
scatter means that adopting optical/UV continuum lu- 
minosity as a proxy of the size of the broad-line region 
introduces an uncertainty of ~ 0.10-0.15 dex in the Mbh 
estimate, with the larger r.m.s. for the line luminosities. 
The velocity scale that best matches the crune of H/3 is 
the line width of Mgll, which gives an r.m.s. scatter of 
0.17 dex. This is better than the CTiino of Ha (0.20 dex) 
and the FWHM of H/3 (0.23 dex) - as expected from the 
large distribution of FWHM/crnne ratios - which are in 
turn slightly better than the other indicators. The worst 
match is obtained for line luminosities with FWHM, with 
an r.m.s. scatter of ^ 0.24— 0.33 dex. From this study, we 
infer a lower limit to the relative accuracy of the various 
indicators of order 0.1-0.2 dex, depending on the choice 
of estimators. If the relationships presented here were to 
be extended beyond the range of Mbh considered, it is 
likely that the scatter will increase, as suggested by the 
slope of the points in the corresponding panels. Other- 
wise, the slope /3 of the size luminosity relation will have 
to be fitted independently. 

6. SUMMARY 

In this paper we have used Keck and SDSS spectra 
of nineteen Seyferts at z = 0.36, to perform a compre- 
hensive study of "virial" black hole mass estimators for 
broad line AGNs. The main results can be summarized 
as follows: 

1. We have fit Gauss-Hermite series to the data in 
order to measure the FWHM and criino of Mgll, 
H/9 , and Ha , as well as Ha and H/3 luminosities 
and continuum luminosities at 3000A and 5100A. 
Measurement errors are approximately 0.02 dex on 
the Mgll and H/3 hne widths, 0.04-0.05 on the Ha 
line widths, 0.01 dex on the continuum luminosity. 
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Fig. 10. — Comparison of Mbh estimates according to the new cross-calibrated formulae as discussed in Section |5] A relation of the 
form log Mbh= <^ + Slogviooo + /^logL is assumed, where viooo is a velocity estimator in units of 1000 km s~^, and L is a luminosity 
estimator in units of 10* erg s~^ or 10*^ erg , respectively for continuum or line luminosity. The slope /3 is fixed at 0.47 for L300O1 
0.518 for Lsioo.n, 0.69 for Lsioo.ti 0.56 for Ljj/Si a-nd 0.55 for Lhq as taken from the literature. The best fit coefficients a are given in 
Table [3] together with the r.m.s. scatter of the comparison. The typical measurement error bar is 0.05 dex. 



0.01 dex on the H/3 luminosity, and 0.06 dex on the 
Ha luminosity. 

2. We have compared twelve formulae taken from the 
literature, showing that Mbh estimates can differ 
systematically by as much as 0.38 ± 0.05 dex (or 
0.13 ± 0.05 dex, if the same virial coefficient is 
adopted). Such differences should be taken into 
account when comparing data obtained with dif- 
ferent methods. 

3. We have cross-calibrated a set of 30 empirical 
recipes based on all combinations of the veloc- 
ity and luminosity indicators corresponding to the 
Mg II, H/3, and Ha broad lines. Taking the masses 
measured by Treu et al. (2007) as our fiducial 
black hole masses, we find that: the absolute scale 
of the different indicators is calibrated to within 
~0.05 dex; the best agreement is found when using 



the line dispersion of H/? as a velocity estimator, 
with the residual 0.1 dex r.m.s. scatter resulting 
from the various continuum luminosity estimators; 
adopting the line dispersion of Mg II raises the scat- 
ter to 0.2 dex; for the other estimators the intrinsic 
scatter is in the range 0.2-0.38 dex. This implies 
a lower limit of 0.1-0.2 dex on the validity of each 
estimator for each individual case. 

The newly calibrated recipes should be useful to reduce 
the sources of systematic uncertainties when comparing 
different studies. 

This work is based on data obtained with the Hub- 
ble Space Telescope - obtained at the Space Telescope 
Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. These observations are as- 
sociated with program 10216 and with the 10m W.M. 
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TABLE 1 
Sample properties 



Name RA (J2000) DEC (J2000) z i' 



(1) 




(2) 






(3) 






(4) 


(5) 


SOI 


15 


39 


16. 


,23 


+03 


23 


22 


.06 


0. 


,3592 


18.74 


S02 


16 


11 


11, 


,67 


+51 


31 


31 


.12 


0. 


,3544 


18.94 


SOS 


17 


32 


03, 


,11 


+61 


17 


51 


.96 


0. 


,3588 


18.20 


S04 


21 


02 


11, 


,51 


-06 


46 


45, 


.03 


0. 


,3578 


18.41 


S05 


21 


04 


51, 


,85 


-07 


12 


09, 


.45 


0. 


,3530 


18.35 


S06 


21 


20 


34, 


,19 


-06 


41 


22, 


.24 


0, 


,3684 


18.41 


S07 


23 


09 


46, 


,14 


+00 


00 


48, 


.91 


0, 


,3518 


18.11 


SOS 


23 


59 


53, 


,44 


-09 


36 


55, 


.53 


0, 


,3585 


18.43 


S09 


00 


59 


16, 


,11 


+15 


38 


16, 


.08 


0, 


,3542 


18.16 


SIC 


01 


01 


12, 


,07 


-09 


45 


00, 


.76 


0, 


,3506 


17.92 


Sll 


01 


07 


15, 


,97 


-08 


34 


29, 


.40 


0, 


,3557 


18.34 


S12 


02 


13 


40. 


.60 


+13 


47 


56, 


.06 


0. 


.3575 


18.12 


S21 


11 


05 


56. 


.18 


+03 


12 


43, 


.26 


0. 


.3534 


17.21 


S23 


14 


00 


16. 


.66 


-01 


08 


22, 


.19 


0. 


.3510 


18.08 


S24 


14 


00 


34. 


.71 


+00 


47 


33 


.48 


0. 


.3615 


18.21 


S26 


15 


29 


22. 


.26 


+59 


28 


54, 


.56 


0. 


.3691 


18.88 


S27 


15 


36 


51. 


.28 


+54 


14 


42, 


.71 


0. 


.3667 


18.80 


S28 


16 


11 


56. 


.30 


+45 


16 


11, 


.04 


0. 


.3660 


18.59 


S29 


21 


58 


41. 


.93 


-01 


15 


00 


.33 


0. 


.3576 


18.77 



Note. — Col. (1): Target ID. Col. (2): RA. Col. 
(3): DEC. Col. (4): Redshift from SDSS DR4. Col. (5): 
Extinction corrected i' AB magnitude from SDSS photom- 
etry. 



TABLE 2 
Measured Properties 



Name 


CMgll 






FWHMMgii 


FWHMh/3 


FWHMhc 


-i'SOOO 










/h/3, [O/I/] 


/h/3, nt 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


SOI 


2260. 


2133. 


1847. 


4418. 


4755. 


3420. 


1.91 


0.74 


1.66 


2.29 


7.24 


0.1 


0.05 


S02 


2914. 


1928. 


2113. 


4000. 


5188. 


3442. 


2.14 


0.36 


1.51 


3.09 


21.38 


0.17 


0.10 


S03 


2147. 


1745. 


1698. 


3345. 


2945. 


2634. 


3.98 


1.69 


2.82 


3.80 


14.45 


0.13 


0.04 


S04 


2253. 


2392. 


983. 


3636. 


3100. 


2617. 


1.86 


1.42 


2.24 


1.35 


6.92 


0.08 


0.06 


S05 


3784. 


3297. 


2686. 


6315. 


5220. 


3751. 


3.39 


2.04 


2.75 


4.37 


16.60 


0.11 


0.05 


S06 


2353. 


1664. 


1382. 


4023. 


4625. 


4306. 


2.75 


0.54 


2.69 


2.00 


7.94 


0.1 


0.13 


S07 


3297. 


2500. 


2531. 


5561. 


4815. 


4326. 


3.72 


2.26 


3.02 


4.90 


15.14 


0.1 


0.06 


S08 


2365. 


1538. 


1015. 


3380. 


3372. 


3017. 


2.29 


1.25 


2.57 


1.12 


4.47 


0.1 


0.08 


S09 


2542. 


2013. 


1462. 


3824. 


2865. 


2710. 


3.16 


0.78 


3.09 


3.80 


15.49 


0.11 


0.04 


SIO 


2897. 


1690. 


2056. 


4532. 


4410. 


3498. 


7.08 


1.11 


3.80 


4.79 


17.78 


0.11 


0.02 


Sll 


2858. 


1590. 


1435. 


4291. 


2733. 


2569. 


3.24 


0.88 


2.45 


2.75 


10.72 


0.08 


0.01 


S12 


2672. 


3213. 


3232. 


4132. 


9005. 


7163. 


4.37 


1.05 


3.24 


5.01 


23.99 


0.06 


0.02 


S21 


3450. 


3172. 


3208. 


6582. 


7681. 


7493. 


2.69 


2.30 


7.24 


9.12 


63.10 


0.07 


0.03 


S23 


3941. 


3196. 


2688. 


7378. 


9700. 


6870. 


3.09 


1.20 


3.55 


3.47 


14.13 


0.1 


0.05 


S24 


3480. 


2886. 


2667. 


6628. 


7864. 


4468. 


2.82 


0.44 


2.95 


3.47 


12.59 


0.1 


0.04 


S26 


3370. 


1862. 


1657. 


6265. 


5451. 


4440. 


1.78 


0.50 


1.58 


2.69 


6.46 


0.07 


0.12 


S27 


2699. 


1609. 


1157. 


3766. 


2567. 


1832. 


2.19 


0.92 


1.82 


2.45 


8.32 


0.17 


0.05 


S28 


3926. 


2313. 


2221. 


8436. 


5116. 


5412. 


2.45 


0.76 


2.29 


1.91 


6.61 


0.09 


0.03 


S29 


2556. 


1744. 


1520. 


4003. 


3190. 


2216. 


2.09 


0.59 


1.74 


2.04 


9.12 


0.1 


0.05 



Note. — Col. (1): Target ID. Col. (2): (Tune of Mgll (km s~^) measured on the line model fit to Keck data. The average error is 0.017 
dex. Col. (3): (Tune of H^S (km s~^) measured on the line model fit to Keck data. The average error is 0.017 dex. Col. (4): uiine of Ha 
(km s~^) measured on the line model fit to SDSS data. The average error is 0.051 dex. Col. (5): FWHM of Mgll (km s~^) measured on 
the line model fit to Keck data. The average error is 0.017 dex. Col. (6): FWHM of H|8 (km s~^) measured on the line model fit to Keck 
data. The average error is 0.017 dex. Col. (7): FWHM of Ha (km s'^) measured on the line model fit to SDSS data. The average error is 
0.040 dex. Col. (8): Rest frame luminosity at 3000 A in 10** erg s^^. The average error is 0.014 dex. The actual error will be dominated 
by variability of order 10%, as for all luminosities listed in this Table. Col. (9): Rest frame nuclear luminosity at 5100 A in 10** erg s"-', 
from Treu et al. (2007). The average error is 0.08 dex. Col. (10): Rest frame total luminosity at 5100 A in 10** erg s^^. The average 
error is 0.014 dex. Col. (11): Rest frame H/3 line luminosity in 10*^ erg s~^. The average error is 0.011 dex. Col. (12); Rest frame Ha 
line luminosity in 10*^ erg s^^. The average error is 0.062 dex. Col. (13): Narrow component of H^ to [OIIIJA5007 fiux ratio. Col. (14): 
Fraction of flux of the H/3 line in the narrow component. 
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TABLE 3 
Mbh Estimator Factors 



^3000 


LslOO.n 


Lsioo.t 


Ljll3 




(0.47) 


(0.518) 


(0.69) 


(0.56) 


(0.55) 



""Mgll 
0"H/3 

FWHMMgii 

FWHMh^ 

FWHMhc 



7.207± 0.052, 
7.458± 0.027, 
7.588± 0.061, 
6.767± 0.055, 
6.803± 0.069, 
6.986± 0.064, 



0.219 
0.112 
0.260 
0.233 
0.292 
0.270 



7.429± 0.039, 
7.680± 0.021, 
7.810± 0.048, 
6.990± 0.045, 
7.026± 0.056, 
7.209± 0.055, 



0.166 
0.090 
0.205 
0.191 
0.235 
0.233 



7.133± 0.045, 
7.383± 0.028, 
7.514± 0.060, 
6.693± 0.053, 
6.729± 0.071, 
6.912± 0.069, 



0.191 7.150± 0.053, 0.226 

0.119 7.401± 0.038, 0.163 

0.253 7.532± 0.074, 0.313 

0.224 6.711± 0.059, 0.251 

0.300 6.747± 0.077, 0.327 

0.293 6.930± 0.071, 0.303 



6.824± 0.051, 
7.074± 0.042, 
7.205± 0.075, 
6.384± 0.056, 
6.420± 0.079, 
6.603± 0.073, 



0.218 
0.178 
0.319 
0.236 
0.335 
0.311 



Note. — Normalization constants for Mbh- Entries are oit error, r.m.s. scatter of log Mbh vs. log Mt (Fig. IIOL The a 
coefficients are determined using the general formula, log Mbh= o + 2 log viooo + /3 log Lj where viooo is the velocity estimator 
in units of 1000 km s~^, and L is the luminosity estimator, which is divided by 10** erg s~^ for continuum luminosity 
measurements and by lO*'^ erg s~^ for line luminosity measurements. The values used for /3 are listed below the luminosities. 



